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Abstract 



.B-meson light-cone distribution amplitudes are discussed in these lectures in the 
framework of HQET. The evolution equation for the leading- twist distribution am- 
CN | plitude is derived in one- loop approximation. QCD sum rules for distribution ampli- 

^ ■ tudes are discussed. 

(N 
(N , 

^§ ■ 1 Introduction 

in 

Many exclusive 5-decay amplitudes in the framework of SCET contain 5-meson distri- 
^3 ' bution amplitudes [T]. The amplitude of the decay B — > at large photon energies is 
given, up to power corrections, by a convolution of a hard part (perturbatively calculable) 



and the B -meson distribution amplitude |2J. Amplitudes of some decays, e.g., B — > ixlv 
at large pion energies, contain both factorizable and non-factorizable contributions [31 H] . 
Factorizable parts of decay amplitudes contain light-cone distribution amplitudes of the 
initial 5-meson and final hadron(s). They describe large-distance (soft) structure of these 
hadrons, and cannot be calculated in perturbation theory. The theory of hadronic distri- 
bution amplitudes in QCD is reviewed in [Hj. 

Quark-antiquark distribution amplitudes of I?-meson in HQET were introduced and 
investigated in [Hj. They are defined as Fourier transforms of matrix elements of some 
gauge-invariant bilocal operators between 5-meson and vacuum. Renormalization of these 
operators was calculated in one-loop approximation. However, an unusual term 1/e 2 in 
the one-loop renormalization constants was erroneously omitted in The correct evolu- 
tion equation for the leading-twist distribution amplitude was derived in [7j at one loop. 
The evolution kernel contains, in addition to terms obtained earlier [U], an unusual term 
log(u> / fi)S(uj — uj'). The method of solution of the evolution equation is also discussed in |7 
in detail. 

Quark-antiquark-gluon distribution amplitudes of 5-meson and their relations to quark- 
antiquark ones (based on equations of motion) are discussed in jH]. 

Sum rules for the quark-antiquark distribution amplitudes were obtained in [0] . A 
simple model of these distribution amplitudes at a low normalization scale (of order of 
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hadronic scale) was proposed. Radiative corrections to the perturbative term and the 
quark-condensate term were later calculated 0- 

In these lectures, we first briefly discuss what is HQET 1 (Sect.EJ). A much more detailed 
presentation can be found in the textbooks [HI Ej. After a short discussion of fs (Sect El), 
quark-ant iquark distribution amplitudes are introduced (Sect.lH). Quark-antiquark-gluon 
distribution amplitudes and their relations to two-particle ones are discussed in Sect. El 
Sect. |B]is the central (and longest) one. Here renormalization of light-cone bilocal quark 
operators in HQET is considered in one-loop approximation. A detailed derivation of the 
evolution kernel is presented, based on the methods of [Ej. Finally, sum rules for the 
distribution amplitudes are briefly discussed in Sect. 



2 Heavy Quark Effective Theory 

Let's consider a heavy antiquark with mass m and momentum 

p = mv + k (2-1) 

in QCD. Here v is a fixed 4- velocity (v 2 = 1), and the residual momentum k, as well as 
momenta of all light quarks and gluons, are supposed to be small (compared to m). The 
propagator of this heavy antiquark can be written as 

m — mi — k 1 — i 1 _ / 1 \ , . 

+ O — . (2.2) 



(mv + k) 2 — m? + iO 2 k-v + iO \m 
The leading term here is the HQET propagator. This can be graphically presented as 2 

* — t ' = g ■ + O (-) . (2.3) 



in 



mv + k mv + k 

In a vertex sandwiched between such propagators, we may substitute 

H/H^^H, (2 .4) 

2 2 2 v ; 2 K J 

Projector can also be inserted near external legs, so that all QCD vertices can be replaced 
by the HQET ones, 

IX a 



% = ig t a (-v») , (2.5) 
up to O(kfm) corrections. 

1 Unlike most texts on HQET, we consider not a heavy quark but a heavy antiquark. Of course, this 
makes no difference, but the active participant in the distribution amplitudes is the light quark (and, 
possibly, a gluon etc.), and this choice makes notation more natural. 

2 The arrow here is somewhat misleading: there is a particle called antiquark and propagating from left 
to right; it has no antiparticle. 
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These Feynman rules can be obtained from the Lagrangian 



L = Q v iv ■ DQ V + (light fields) (2.6) 
The antiquark field is our main field here, it satisfies 

Qvi> = -Qv\ (2.7) 
Q v is the conjugate field for Q v . Here the covariant derivatives are 

D„q = (^ - iA„) q, qD^ = q (t ^ + , A„ = g A a 0fl t a . (2.8) 



QCD tree diagrams are reproduced by HQET up to 0(ki/m) corrections (they can also be 
reproduced, if we add the appropriate 1/m terms into the Lagrangian). 

The heavy quark chromomagnetic moment is ~ 1/m by dimensionality. At the leading 
order in 1/m, the heavy-quark spin does not interact with the gluon field. Therefore, it 
may be rotated at will, without changing the physics (heavy-quark spin symmetry). It 
may even be switched off (superflavour symmetry). We shall work with the spinless heavy 
antiquark, 

L = Q* v iv ■ DQ V + (light fields) , (2.9) 

during most of these lectures, because this greatly simplifies reasoning and calculations. 
Here again Q* v is the main (scalar) field, and Q v is its conjugate. 

So, tree QCD diagrams, expanded in ki/m to some order, are reproduced by the cor- 
responding HQET diagrams. But what about loops? Here things are not so simple [T3] . 
Let's consider, for example, the heavy-light two-point diagram (Fig. with p = mv + k, 
where the residual momentum k is small. By choosing v along p we can always ensure 
k = ujv. Let's consider the integral 



-im 2 f d d l 



K d/2 J [m 2 - (mv + k + I) 2 ] 2 (-Z 2 ) 



(2.10) 



which has neither ultraviolet (UV) nor infrared (IR) divergences. There are two regions of 
the loop momentum I in this integral: 

• Hard region I ~ m. Expanding the integrand in C m, /, we have 

1 +4 (m + ^ } " + -.. (2.11) 



\m 2 



(mv + I) 2 } 2 (-1 2 ) [m 2 - (mv + Iff (-1 2 ) 



Soft region I ~ uj. Expanding the integrand in k, I <C m, we have 

1 (k + l) 2 . 

+ 2 r - ,} r - ' l3/ - + ••• (2.12) 



-2m(k + I) ■ v] 2 H 2 ) [-2m(k + I) ■ v} A (-/ 2 ) 
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p p + 1 p 
Figure 1: Heavy-light two-point diagram 



The contribution of the hard region is 



I h = m 
h 



-2e 



UJ 



M(2, 1) + 2— [M(3, 0) - M(2, 1) + 2M(3, 1)] 
m 



P *h 



l,m/l m 
-— +log-+ _-21og— - 
r(l + e) 2e: p \e p 

where the on-shell massive two-point integrals (Fig. |2J 

d d l 



1 I - + ■ 

m 



m d/2 m d - 2{ni+n2) M(n u n 2 ) 



arc 



D 1 = m 2 -(l + mv) 2 - iO , £> 2 = -Z 2 - *0 
r(d - ni - 2n2)r(-d/2 + m + n 2 ) 



r(ni)r(d - ni - 7i 2 ) 




M(Tii,Ti 2 ) 

This contribution is IR divergent. 



mv I + mv mv 
Figure 2: On-shell massive two-point integrals 
The contribution of the soft region is 
I* = 



UJ 



-2co)- 2£ 1(2, 1) + — [1(3, 1) - 2/(2, 1)] + 



P 2£ l 



1 -2u 
r(l + e) ~ 2e ~ ° g ~p~ 



1 nl -2cu 1\ 

__21og - + 

e p 2 J m 



where the HQET two-point integrals (Fig. ^ u = k ■ v 

d d l 



d/2, 



2uj) d - ni - 2n2 I(n 1 ,n 2 ) 



are 



D x = -2(1 +p) -v-iO, D 2 = -l 2 -i0 
, T(-d + n 1 + 2n 2 )T(d/2-n 2 ) 




k l + k k 



Figure 3: HQET two-point integrals 

This contribution is UV divergent. 
The complete result is finite: 

/= -, og =^ + (21og^-iW... (2.19) 

This expansion can be easily continued if desired. 

What about higher loops? Let's consider, for example, the two-loop two-point diagram 
with a small external residual momentum. There are several regions of the loop momenta 
in this diagram (Fig. HJ. In each of them, some momenta are hard (~ m), some are 
soft (~ lu). As we have already seen, a heavy-quark line with a soft residual momentum 
becomes an HQET line. Soft massless lines are shown by dashed lines in the figure. 




Figure 4: Regions in the two-loop diagram 

For example, in the second diagram in Fig. the left loop is hard. It contains a single 
scale m. All external momenta of this loop are soft, including those of the lines belonging 
to the other loop. We can expand the integrand in these small momenta; after taking the 
loop integral, we obtain a polynomial in these momenta. Now from the point of view of 
the soft loop (large distances) this hard loop is just a local vertex. This loop contains a 
single scale u, and we can calculate it, obtaining a non-analytical function of the external 
residual energy ui. 

In a general multiloop diagram, hard lines must always form loops, so that momentum 
conservation can hold after neglecting all soft momenta. There can be several disconnected 
hard parts; each of them must contain at least one heavy line (a subdiagram consisting of 
only light lines and having soft external momenta has no reason to be hard). From the 
point of view of the soft part, hard parts are just local vertices. There may appear a soft 
subdiagram connected to the rest of the diagram only at such a vertex. Such a subdiagram 
is scaleless and hence vanish. For example, we could consider one more region in Fig. 
when the only soft line is the middle light line, all the rest are hard. But this soft line 
forms a loop containing one local vertex (the integrated hard loop), and hence it vanishes. 
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In the usual HQET formalism, Lagrangian contains local operators multiplied by match- 
ing coefficients. QCD operators are also expanded in HQET operators with matching coef- 
ficients. These matching coefficients are the only quantities in the theory which depend on 
the hard scale m. Diagrammatically, they come from hard loops in QCD diagrams. Local 
operators produce vertices polynomial in their external momenta. They appear in HQET 
diagrams, which contain only the soft scale u. These HQET diagrams are soft parts of 
QCD diagrams. 

3 B-meson decay constant 

During these lectures, we shall mostly live in a world with a heavy antiquark having 
j p = + . Physics in the real world is the same, up to 1/m corrections. We shall work in 
the v rest frame. 

The ground-state S"-wave Qq meson has the quantum numbers j p = | + . There are 2 
P-wave excited mesons with j p = | and | . The heavy-light quark current 

3 = Qti (3-1) 
has no definite parity; the currents with parity P = ±1 are 

Jp = g 3 ' ( ) 

They have the quantum numbers of S-wave | + mesons and P-wave | mesons. The 
ground-state meson M has a Dirac wave function u which satisfies 7°w = u and is normal- 
ized by uu — 1. The matrix element of j from M to vacuum is 

<0\j\M>=Fu, (3.3) 

where the one-meson states are normalized by the non-relativistic condition 

<M,p'\M,p> = (2Ttf5(p'-p) . (3.4) 

The correlator of the heavy-light currents (Fig. EJ), 

i<Tj(x)j(0)> = 8{x)U{x°) , (3.5) 

has the structure 

U(x°) = A + Bf). (3.6) 
Therefore, the correlators of the currents with definite parity are 

i<Tjp(x)j P (0)> = S(x)U P (x°) , Up = A + PB = ^Tr(l + P 7 °)n. (3.7) 

The spectral density of, say, the correlator with P = +1 is 

p+(e) = F 2 5(£-A) + --- (3.8) 
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x 



Figure 5: Correlator of heavy-light currents 



where A is the residual energy of the ground-state meson, and the dots mean contribution 
of excited states. 

Now we shall return to the real world with j p = | heavy antiquark for a while, but 
still with m = oo. The 5*-wave ground-state meson turns into a degenerate doublet with 
j p = l - . The two P-wave mesons turn into two degenerate doublets: with j p = + , 
1 + and with j p = 1 + , 2 + . All heavy-light currents Q v Fq reduce, due to (|2.7jl . to 4 ones 
with 

T = 75 , 7 and T = 1 , 775 . (3.9) 

The first 2 currents, with F anticommuting with 7 , have the quantum numbers of the 
ground-state _ , 1~ mesons; the second 2 currents, with F commuting with 7 , have the 
quantum numbers of the P-wave + , 1 + mesons. 
The correlators (Fig. EJ) are 

i<Tj 2 (x)jt(0)> = 6{x)U 12 {x°) , n 12 = Trrx^^ran, (3.10) 

or 

-1-7° 

n 12 = n P Trr 1 ^^r 2 , (3.11) 

where P = +1 for F anticommuting with 7 and —1 for commuting. For T = 75 and 7*, 
the correlators are 2H+ and 2Fl + 5 l i . Their spectral densities are F^5(e — A) + • • • and 
Fz*5(e — A)8 ij -\ , where 

<0\Q vl5 q\B> = F B , <0|Q,7g|P*> = F B *e , (3.12) 

and dots mean contribution of higher states. Therefore, 

F B = F B *=V2F. (3.13) 

The usual definition of the decay constants is 

<0\Q v Yl 5 q\B> r = if B p» , <0\Q v Yq\B*> r = imf^e? , (3.14) 

where the relativistic normalization 

r <B,p'\B,p> r = (2ir) 3 2p°5(p - p) (3.15) 
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of single-meson states is used (this normalization becomes meaningless in the limit m — > oo, 
and therefore cannot be used in HQET). Comparing ()3.14|) for B or B* at rest with (|3.12jh 
we arrive at 

2F 

f B = f B . = —, (3.16) 



up to 1/m corrections. 

Let's consider the correlator ()3.11|) with P = +1. Equating the ground-state contri- 
butions to the spectral densities of the left-hand side and the right-hand one, we obtain 



Therefore, 



1 + -Y - 

<0\j 2 \MxM\j+\0> = F 2 Tr r 2 ^-— r\ . (3.17) 



<0\Q v Tq\M> = ^Trr^^r M , (3.18) 



where the matrix 

\if for B* 

is defined up to a phase factor. We can re-write this result for the relativistic normalization 
of the meson state, in covariant notation: 

<0|Q„rg|M> r = ^/mF Tr YM , (3.20) 

where 

M= 1 -±^h 15 iOTB (3.21) 

2 for B* y ' 

(of course, one can re-define the phases of \M> and hence of M). 



4 Quark— ant iquark distribution amplitudes 



After safely returning to the ideal world with a + heavy antiquark, we want to invent an 
operator which probes more details of the structure of I?-meson than the local current (J3.1j) . 
To this end, we consider a bilocal gauge-invariant operator 



6(t) = Q* v (0)[0,z]q(z) 



where 



and 



[x, y] = P exp 



z = , t = v ■ z . 



(4.1) 

(4.2) 
(4.3) 
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Its matrix element from the ground-state meson to vacuum has 2 Dirac structures: 



<0\6(t)\M> = F 



~ ^ , <P-(t) - <P+(t) 4 

¥>+W + 2t f 



u. (4.4) 



because fu — u. 

In what follows, we shall often use light-front components of vectors. Let's introduce 
(in the v rest frame) two light-like vectors 



n£ = (l,=Fl,0), 

n 2 , = n 2 _ = , n + ■ n_ = 2 



(4.5) 



Light-front components of any vector a are defined as 

a± = a ■ n± = a ± a 1 . (4-6) 

We have 



a ■ b = - (a + 6_ + a_fc+) — a± ■ bj 
2 



(4.7) 



In particular, 

= ^-(n$+n l L) , v+ = v- = l, v ± = 6. (4.8) 
We shall also use light-front components of 7 M : 

7± = 7 • n± = fi± ■ (4-9) 
The definition ()4.4)1 can be re-written as 

<0|O(t)|M> = ±F[<p+{th- +<P-(th+]u. (4.10) 
If we introduce the operators 

± (t)= 7± 0(t), (4.11) 

then 

<0\O±(t)\M> = F^±(*)7±u. (4.12) 
The 5-meson distribution amplitudes are the Fourier transforms of these functions: 

cp ± (u>) = ^- [fi±(t)e wt dt, <p±(t) = f ip±{u)e- iwt du . (4.13) 



2vr „ 

They are normalized by 

£ ± (0) = / v± (uo)du: = 1. (4.14) 
9 



The function <p+(w) is the leading- twist distribution amplitude, and <p-(uS) - the subleading- 
twist one (though there is no good definition of twist in HQET). We can formally introduce 
the operators 

f J (4.15) 
6 ± (t) = / 0±(uj)e~ lu)t duj , 

then 

<0|O±(u;)|M> = Fy? ± (cu)7 ± M . (4.16) 

The distribution amplitudes describe the distribution in the light-front component p + of 
the light-quark momentum in i?-meson. 

The expansion of the operators (j4.11j) in t reads 



6±(t) = 5> 



(n)(-it) n 



n=0 

Of = J ± (uj)uj n du = Ql 1± (iD + ) n q, 

or for matrix elements 

-it) n 



(4.17) 



<p±(t) = J2 <ujU> ±~ 



n< 

n=0 



<u n >± = / ip±(u)u n du}, 
Jo 

<0\O ( ± ) \M> = F<uj n >±-f±u. 



(4.18) 



We can also reconstruct 0±(u) from O^: 



± (u)=l 0^—'^, (4.19) 



(00 



or for matrix elements 

r +tco dn 
¥>±(w) = / <w n >±o;- n - x — (4.20) 

(you can easily check this by substituting ()4.17|) into (|4.19j) : integration in dn yields 
5{\og{u'/u))). 

The first moments <u>± can be found from the equations of motion. The equation of 
motion for the heavy antiquark is 

Q* v D o = 0, (4.21) 

and therefore we obtain 

<0\Q* v D q\M> = <0\d {Q* v q)\M> = -iFku. (4.22) 



10 



The vector part has the structure 

<0\Q* v Dq\M> = aF^u . (4.23) 
The equation of motion of the light quark is 

<0\Q* v l/)q\M> = 0, (4.24) 
and using I/) = -Dq7° — D ■ j we obtain 



<0\Q* v Dq\M> = -FAju . (4.25) 
Finally, we arrive at 

4 _ 2 - 

<cu> + = -A, <u>_ = -A. (4.26) 

3 3 

Let's now consider the second moments. They involve two new (nonperturbative) 
hadronic parameters: 

<0\Q* v E ■ aq\M> = —%F\ 2 E u , <0\QIH ■ aq\M> = —F\ 2 H u , (4.27) 

where 

E = i[D ,D], H = iDxD, a = 7 °7, a = -7757 . (4.28) 

Now we can calculate all matrix elements with 2 derivatives. From (|4.22|) and ([4.250 . using 
the heavy- ant iquark equation of motion ([4.21)1 . we immediately find 

<0\QlD 2 q\M> = -FA 2 u , <0\Q* v D Dq\M> = \fA 2 ^u . (4.29) 

Using the definition ([4.27)1 of X 2 E , we immediately find 

<0\Q* v DD o q\M> = ^F (A 2 + A|) iu . (4.30) 

The second spatial derivatives have the structure 

<0\Q* v D i D J q\M> = F (b6 ij - ^X 2 H e ijk a k j u , (4.31) 

where the second coefficient follows from the definition ([4.270 of \ 2 H . We find h using the 
light-quark equation of motion ([4.240 : 



<0\Q* v D i D j q\M> = F 



(A 2 + A| + X 2 H ) & + % -\\^ k o k 



u. (4.32) 



Finally, we arrive at 



<u 2 > + = 2A 2 + + ~\ 2 H , <u; 2 >_ = h 2 + ±\% . (4.33) 
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What about 5-meson distribution amplitudes in the real world with a \ heavy anti- 



quark? It has 4 distribution amplitudes, as any pseudoscalar meson: 

<0\Q* v (0)[0,z] l5 q(z)\B> r = -tf B m<p P} 
<0\Q* v (0)[0,z}^ l5 q(z)\B> r = f B [i0 AlP » - m^ A2 z' 
<0\Q:(0)[0,z]a^ l5 q(z)\B> r = ifsMP^" 

Similarly to ()3.2()|) . we have 

<0\Q* v (0)[0,z}Tq{z)\M> r = FTtT <p + + 



(4.34) 



2t 



M 



(4.35) 



Therefore these 4 QCD distribution amplitudes can be expressed via 2 HQET ones: 

<p+(t) + <p-(t) „ „ 

VP = n ' V^Al = <P+{t) , 



(f A2 = (fi T = -- 



i <p+(t) - <p-(t) 



(4.36) 



t 



The QCD distribution amplitudes are usually considered as functions of the longitudinal 
momentum fraction x = u/m. It is usually assumed that (p A i(x) ~ % at x — > 0, and 
<pp(x) — > const. Therefore, we shall assume that <p+(lj) ~ to at uj — ► 0, and (p-(u) — ► const. 
The QCD distribution amplitudes are normalized as 



£p(0) = <p A1 (0) 



(p A2 (0) = ip T (0) 



A 
3 



(4.37) 



(to derive the last formula, we used the t expansion ()4.18|) and the first moments (|4.26j) ). 

The vector meson B* is described by 6 distribution amplitudes in QCD. All 10 distribu- 
tion amplitudes (4 for B plus 6 for B*) are expressed via 2 HQET distribution amplitudes 
(p±(u). This is a consequence of the heavy-quark spin symmetry. We don't present formu- 
lae for B* here; they can be found in |B]. 



5 Quark— ant iquark— gluon distribution amplitudes 

Now we shall discuss relations of quark-antiquark distribution amplitudes and quark- 
ant iquark-gluon ones following from the equations of motion. In order to apply them, 
we need to differentiate with respect to the coordinates of the light quark and the heavy 
antiquark separately. Therefore, we go slightly off the light cone. The generalization 
of (|4.4j) to an arbitrary z 2 is 



<0\Q;(0)q(z)\M> = F 



z ) + 2t * 



u . 



(5.1) 



where t = v ■ z, and the fixed-point gauge x^A^{x) = is used to simplify notation. 
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In order to apply the light-quark equation of motion (|4.24J1 . we apply 

1 dz» 

to this definition. The differentiation of the right-hand side is straightforward. In the 
left-hand side we write 

<0\Q* v (0)Y(d, - iA^z) + iA^z))q(z)\M> . 
The first two terms yield zero. In the fixed-point gauge 

A li (z)= [ G u ^{uz)uz u du, (5.2) 
Jo 

where G^ = gG a ^t a - 

The matrix element of the operator containing G^ contains 4 structures: 



<0\Q* v (0)[0,uz)iG 1/fl (uz)z' / [uz, z]q(z)\M> = 



u . 



(5.3) 



and is parametrized by 4 quark-antiquark-gluon distribution amplitudes (note that this 
matrix elements vanishes when multiplied by z M ). 

Equating the coefficients of u and jku, we obtain two consequences of the light-quark 
equation of motion: 

~ _ ~ r 1 

+ f__P± = 2t (jpA - ^v)u du , (5.4) 
* Jo 

(p' + - #_ + + 4 ^ = 2t J o ( 2 Vv + 4>a + i>x)u du . (5.5) 

Similarly, to use the heavy-antiquark equation of motion we apply 

Q* v (0)v^q(z) = v^(Ql(0)q(z)) - v»Q* vX Q)(d » + iA^O) - iA^))q{z) . 

The first two terms in the last operator yield zero. Now we move the heavy antiquark; 
the center of the fixed-point gauge must be constant, and we have to use the gauge (x — 
z^A^x) = 0: 

= - / G^{uz){\ - u)z v du . (5.6) 
Jo 

We arrive at two consequences of the heavy-antiquark equation of motion: 

<P+ + ^ ~ ^ + + *A^+ + 2t^± = ~tj {4>A + 4>x)(l-u)du, (5.7) 
Cp + — (p- - „ . f 90- dip 



0--0+ + z f + iA(<p--<p + ) + 2t[^--^j =2tJ (ij> A +ifr)(l-u)du. 

(5A 
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The functions 

dfi±(t,z 2 



,2=0 

are some new (non-leading) quark-antiquark distribution amplitudes. We are not inter- 
ested in them here. There are 2 equations involving only our familiar distribution ampli- 
tudes <p±(t, 0), namely, (|5.4|) and a combination of ()5.5|) and (|5.7jl : 

+ - — — — = 2t (ip A - ipv)udu , 

ri . . (5 " 9) 

+ </?'_ + 2iAy> + = -2t / (tp A + ip x + 2ip v u)du . 

Jo 

The quark-antiquark-gluon distribution amplitudes in the momentum space are defined 

as 

4>i(t,u) = J \pi{uj,Oe~ i{u}Hu)t dujd^. (5.10) 
Performing Fourier transform of the equations (J5.9)) . we obtain 

w -iT + ^ w ) = I ( w )' (5.11) 
(c<j — 2A)(/) + (u;) + c<j<y9_(u;) = J(c<->) , 

where 



1 dp L jdt^P'®-^^ 



d f u f°° d£ C w f°° d£ 8 

J(uj) = -2— dp fWM) + fe(M)]-4 dp ~<MM)- 
duJ Jo Ju-p € Jo Ju-p i o£ 

If we insist on having ip+(0) = 0, then the condition 

f°° d£ 

j(o) = -2/ -f [Mo,Z) + Mo,Z)] = Q ( 5 - 12 ) 

Jo s 

must be satisfied. It should follow from vanishing of ipA,x(^,0 at u; — > (behaviour of 
the three-particle distribution amplitudes at the boundaries is discussed in [T4~]). 

If we new iI>a,v,x{w, £)> we could solve the equations (|5.11j) for </?±(a;). This is not very 
realistic, because we don't know them. The solution is a sum of two terms: 

<p ± {u)=<p<r^(u,)+<p%\u,), (5.13) 

where (p^± W \uj) is the solution of (|5.1H1 in the case if all quark-antiquark-gluon distri- 
bution amplitudes vanish (it is called the Wandzura-Wilczek part of the solution), and 
(p±\u) is induced by the gluonic terms. The Wandzura-Wilczek part is 

*>r>) = 5^(2A--), *£^M = ^»(2A-«) (5.M) 
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(Fig. 1^1); it satisfies the normalization conditions (j4.14|) . Note that 2A is the maximum 
value of p + of the light quark if we assume that 5-meson (having the residual energy A) 
consists of the on-shell heavy antiquark (always having zero residual energy) and the on- 
shell light quark. The gluon-induced part is given by some explicit integrals of ipA,v,x(&, 0; 
we don't present these long expressions here, they can be found in [Hj. 




Figure 6: The Wandzura-Wilczek parts of y?+(co>) (solid line) and (p-(ou) (dashed line) 
The moments 1)4.18)1 also consist of two contributions. The Wandzura-Wilczek parts 



are 



(W) _ 2(2A) W n (w) _ 2(2A) n 



+ n + 2 ' _ (n + l)(n + 2) v ' 

The gluon-induced part does not contribute to the zeroth moments (normalization) and the 
first ones; its contribution to the second moments (J4.33)) is expressed via the normalizations 
of the quark-antiquark-gluon distribution amplitudes: 

J Mu,0dud£=l>?H, (5-16) 



The contributions to the higher moments are expressed via the moments of ?Pa,v,x(w, £); 
the explicit formulae can be found in [S]. 



6 Evolution 

Until now, we neglected renormalization of the considered operators. In this section we 
shall discuss renormalization of the leading-twist 5-meson distribution amplitude. 
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Let's summarize what we know about the bare operators of interest. There are 3 
families of such operators, or 3 "representations" (t, u, n): 



O + (t) = Q* v (0) 1+ [0,z]q(z), 
O + (u)=Q* v (0) 1+ 8(iD + -u;)q(0), 



(6.1) 



O 



g;(o) 7+ (^+) n g(o). 



They can be converted into each other: 

r 00 i 

6 + (t) = / + {uo)e- iujt duo = Y,°+ } - 



-it) n 



Oju) 



+ {t)e 



IU!1 



dt 
2^ 



O 



(n) 
+ 



dt 



o+(t) 



n=0 

27T? 

oo 

+ (uj)uJ n duj . 



(6.2) 



t=o 



As we shall see, not all of these relations survive renormalization. 

We shall calculate matrix elements of these bare operators, therefore, we need the Feyn- 
man rules for them. If we retain id + in all brackets (iD + ) n in 0+ , we obtain the quark- 
antiquark vertex with shown in Fig. [TJ When transformed to the ^-representation (|6.2j) . 
this gives 5{p + — uj), as expected from the form ()6.1Jl of + (u). 



a 



to 




P 



P+7+ 




P 



8(p+ - w)7+ 



Figure 7: Quark-antiquark vertices 
Now let's retain a single A + in (iD + ) n . After some combinatorics, this gives 

.)"■ (6.3) 



id + + A + r - jr ( n m \ [(id+^A+i {id, 

m=l ^ ' 



This gives the quark-antiquark-gluon vertex of 0+^ with 



E 

m=l 



n 
m 



Um-l^n-m _ (P+ + ^+)" ~ P+ 



(6.4) 



shown in Fig.|HJ Transforming it to the cu-representation ()6.2|) gives two 5-functions (Fig.|HJ). 



Of course, the integral in duj of this vertex vanishes, because 0+ does not interact with 
gluons. 
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n 




ji a 



{P+ + K) ~ P+ .a v 



8{p+ +k+ -a)- 5(p+ - u) ^ +a M 



Figure 8: Quark-antiquark-gluon vertices 

The bare operators + {uS) can be expressed via the renormalized operators + (to; fi) 
(n is the MS renormalization scale), or vice versa: 

(6.5) 
(6.6) 

The renormalization "matrices" Z + (to,to'; fi) and Z^ x (uj, to'; ji) are inverse to each other: 



+ (to) = J Z + (ou,u'; fi)0 + (u'; fi)doo' . 
+ (to; fi) = / Z+ 1 {lo,lo'\ jj)0 + {uj')du' . 



Z + (to, to"; n)Z + (to", to'; fi)du" = 5(uo — to') , 
(u, to"; fi)Z + (to" , to'; n)dto" = S(to — to') . 

The renormalized operators + (to; fi) obey the renormalization-group equation 
dO + (uj;fi) + f r+ ( WjU / ; ^ o+ ( w / ; ^ tfc; / = 0j 



(6.7) 



d log fx 

where the anomalous dimension "matrix" is 



(6.8) 



,, , / / ry — i / // \9Z + (u" ,u'; fl) „ 
r + [to,to ; fi) = / Z + (to, to ; fi) — duo 



d log /! 



dZ^{to,uo";^) 
dlogfi 



(6.9) 



Z + (u> , to'; jj,)duo" . 



This equation can be derived in two ways: either we differentiate (|6.5j) in d log fi and obtain 
(because the bare operator is /i-independent), or we differentiate (|6.6|) . 
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At one loop 



where 



Z + (u, u>'; fj) = 5(u — u') + 2+ (o>, lo'\ fi)a s + ■ 
Z+ l (u, uj'\ /j,) — 5(u — u) 1 ) — z+ (u, uj'\ fi)a 8 + • 
r + (u,u';fi) = r^(u,u';fi)a a H 

as = —- 



(6.10) 



From (|fj.9jl we obtain 



rV = ^ 2ez\'(u,u;n) 

a log fi 



(6.1i; 



The matrix element of 0+(u) between a state with the light quark with momentum p 
and the heavy antiquark with momentum p' (which are off-shell) and vacuum is 



M = <0\O+(u)\q(p), Q:(p')> = Z^Z 1 ' 2 [6(p+ - + Ml + M 2 + M 3 



<0\O + (u; n)\q(p),Qt(j/)> + a s I z$\u, u'; /i)5(p + - u') 1+ du' 



(6.12) 



where M\££ are the contributions of the one-loop diagrams shown in Fig. El the matrix 
element of the renormalized operator is finite at e — > 0, and z+(u,u';fi) contains only 
negative powers of e (in the term with z+ , we may substitute 8(p+ — oj') r i+ instead of 
<0\O + (u'; n)\q(p), Q*(p')>). This allows us to find the renormalization "matrix" with 
one-loop accuracy. 







Figure 9: Diagrams for the matrix element of + (u) between a quark-antiquark state and 
vacuum 

Let's calculate the first diagram (Fig. ITUj) for 

p+ = u', p± = 0, p 2 = p + p- < . 
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It is 



M x = iC F g 2 Q 



d d k S( P+ - uj) - 5(k + - uj) 



(2ir) d p + - k + 

The numerator can be simplified as 

7+^7+ = 2/c+7+ • 



[-(p-k) 2 -i0] [-k 2 -iO] 




(6.13) 



Figure 10: The first diagram 
This diagram can be written in the form 



M x = 2C F - 



fi(u),u') - S(uj - uj') I fi(u",u')du" 



(6.14) 



( 4vr )d/2 

where the term with fi(u, uj') comes from the second 5-function in ([6.130 . This function is 

5(k + — uj) 



7T 



r -/a ( _ p a ) -« r(e)/l(a;jt y ) = _ i _ 
Now we use a parametrization 



UJ 



UJ — UJ 



<fk 



[-(p- A;) 2 -20] [-k 2 -i0] 



1 



-k 2 - iO 



i^+ i0 ^da 



for both denominators, and also write the ^-function as 



dv 

S(k+ -uj) = 25(2(k+ -uj)) = 2 exp [2(fc+ - — . 

J-oo 27T 



We obtain 



^(V)-T(e)/i(w,oO = -2<- 



/d d k da\ da2 — 
Z7T 



x exp [oti(k — p) 2 + 012k 2 + 2%v{k ■ n + — uj 
Shifting the integration momentum as 



k' = k 



aiip — tvn + 



(6.15) 



(6.16) 
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to form a full square, we obtain 

CO 



-2i- 



u' — UJ 
dv 



x , -exp 



The integral in d d k! is calculated using the Wick rotation k = ikEo'- 



da.\ da2 exp 
2i, • ai 



aia 2 o 

— : — P 

Oil + 012 



-UJ — UJ 



d d k'exp [(«i + a 2 )(k' 2 + iO 



J d d ke a ^ + ^ = i J d d k E e- ak * = % 



n y/2 



(6.17) 



The integral in dv gives a 5-function. We obtain an integral in two a-parameters: 

(-p'r^Muj^') 



UJ 



UJ' — UJ 



J da\ da 2 («! + a 2 ) d ^ 2 exp 



ai«2 2 
— T — P 



5 



Oil , 

-UJ — UJ 



oci + a 2 



It is always possible to calculate one integral in a "radial" variable rj in the space of a- 
parameters via T-function. In this case, the most convenient choice of such a variable is 
rj = ai + a 2 - Therefore, we insert 5(o>i + a 2 — r])dr] under the integral sign, and substitute 
ai = rjXi-. 



UJ 



UJ — UJ 

The final result is 



dxi dx 2 S(xi +x 2 - 1)5{xiuj' - u) J d7]7]~ 1+e e- { " p )xiX2V . 
6{uj'-uj) uj^ 



Functions F(u,u') which appear in the evolution kernel should be understood as dis- 
tributions: they are always integrated with smooth test functions f{uj'). The distribution 
[F(uj,uj')] + is defined by 



[F(uj, uj')} ip{u')duj' = / F(uj, uj') {ip{u') - ip{u)) duo' . 



Therefore, formally we can write 



(6.19) 



(6.20) 



The result ()6.14|) . ()6.18|) of the calculation of the diagram in Fig. El can be written via a 
+- distribution as 



F(uj, uj') = [F(uj, uj')} + + 5(uj — uj') / F(u, uo")duo" . 



[fi(u,u')] + +5(u-u') / f\(uj,uj")duj" - / fi(u",u)du" 



(6.2 
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The coefficient of 8{u — uj') here can be calculated by substitution x = uj" /to: 

• i 



x 



-l+2e, 



l-xY l - £ dx- I x l ~ E (l-xY l - e dx 



l-ei 



-Is. 



'1 JO 

Substituting x — > 1/x in the first integral, we have 



[ (x- £ -x 1 ~ £ )(l-x)- 1 ~ £ dx= [ x~ £ {l-x)~ £ dx. 
Jo Jo 



Therefore, the final result for Mi (Fig. fTUj) is 

- uj) u l - £ 



(47r) d / 2 



[uj'-UjV+c (UJ 1 ) 1 -* 



+ 



r 2 (i- g ) 

T(2 - 2e) 



5(cj — a;') 



Now we shall calculate the second diagram (Fig. ITTj) for 

p ■ v = u\ < . 



It is 



M 2 = -iC F gZ 



(6.22) 



(2tt) c 



[-P - iO] [-{p' - k) ■ v - iO] ' 

(6.23) 




Figure 11: The second diagram 
This diagram can be written in the form 



M 2 = 2C F - 



9l 



/ 2 (u; — uj') — 5(uj — uj') / f 2 (uj")du! 



(6.24) 



( 4?r )d/2 J 

where the term with f 2 {uJ — w') comes from the first 5-function in ()6.2H|) . This function is 



^ 2 T{e)f 2 {uj") 



2uj" 



d a k 



[-k 2 - iO] [-(p 7 - Jfe) • u - iO] 



Now we use a-parametrization ()6.15|) for both denominators (it is convenient to multiply 
the linear denominator by 2 before this) and (j6.16|) for the 5-function, and obtain 



-— / d d kda\ da 2 — exp \a\k 2 + 2a 2 (p — k) ■ v + 2iv(k 

UJ J ZTl 



n< — uj 



")]■ 
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Shifting the integration momentum as 



k' = k 



a 2 v — ivn + 
oti 



to form a full square, we obtain 



2i 



OJ" 



— doti da 2 exp 



h 2oj x a 2 

ati 



dv 



x / — exp 



2iv ( oj 



d d k'exp [a^k' 2 + iO)] . 



Taking the integrals in d d k' f)6.17|) and dv, we obtain the integral in two a-parameters: 

T(e)f 2 (oj") = — - I dct\ da 2 exp 
oj J 



h 2uj x a 2 

a,\ 



r|«2 „ 

o C<J 



Now the best choice of the "radial" variable is a±, so we simply substitute a 2 = ot-iy: 

1 

up' 



Finally, 



- / dy5(y-uj") J da ia -{ 1+£ e- y{y - 2uJl)ai . 
f f ^ 9 ^ 



(uj") 1+£ (uj" - 2cji) £ ' 
Now we rewrite M 2 ()6.24|) . ()6.25|) via a +-distribution: 



(6.25) 



M 2 = 2C F - 



rr(e) 7 4 



( 4vr )d/2 

The coefficient of 5(uj — c</) is 



[f 2 (uj - oj')] + +5(oj-oj') ( / f 2 (uj - uj")doj 



(6.26) 
(6.27) 



We introduced tui only to regularize possible infrared problems; we only need the UV 
divergence of this diagram, which does not depend on oj\. Therefore, we may assume 
|a>i| <C uj. The coefficient of S(u — oj') is 



"^-^dw" 



UJ 



-2s 



2e 



and the final result for M 2 (Fig. ITTj) is 



M 2 = 2C F 



9l 



(4vr) rf /- 



■r(e)7H 



6{uj- uj 1 ) 
(uj - uji) l + 2 z 



OJ 



-2c 



2e 



-5(oj - oj') 



(6.28) 



(6.29) 



This one-loop diagram contains al/e 2 UV divergence! It is in the coefficient of 5(uj — uj') 
given by the integral (|6.27p . The function f 2 (oj") has al/e UV divergence in the integration 
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in transverse momentum; the longitudinal integral ()6.27|) is again UV divergent (|6.28|) . 
The operator at the vertex of Fig. ^2 is (after integration in uj) a light-like Wilson line; the 
heavy-antiquark propagator is a time-like Wilson line. The vertex correction to a time-like 
- light-like cusp on a Wilson line is known to have a 1/e 2 UV divergence at one loop [TK] . 
The third diagram (Fig. IT2*|) is 

M _ r 2 [ ddk 5(k+-u) 

3 F9 ° J (2ir) d [-(k-p) 2 -t0}[-k 2 -tO}[-(p' +p-k)-v-tO] ' 1 J 

The numerator can be simplified as 

i+M = k +i+ ■ 

We shall now demonstrate that this diagram is UV-finite and hence does not contribute to 
the renormalization constant. 




k — p 



Figure 12: The third diagram 



M 3 = — AiCpg Q u} r yj r I d kda\da 2 da 3 



Using a-parametrization 1)6.15)1 for the denominators and (J6.16)) for the 5-function, and 
obtain 

dv 
27 

x exp [ax (A; — p) 2 + 02k 2 + 2a 3 (p + p — k) ■ v + 2%v{k ■ n + — uSj\ . 
Shifting the integration momentum as 

ctjjo + a 3 v — iun + 



k' = k- 



we obtain 



— AiCF9Q^1+ J dai da 2 da 3 e~ 



dv 
2^ 



exp 



. a x J + a 3 
liv I UJ 



d d k' 

(2vr) d 



da± da 2 da 3 5 ( — — — — - — to) e A 

Oil + OL2 
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where 



a 3 (a 3 - a 2 uj') + a 2 (a 1 u' + a 3 )(-p 2 ) , . 
A = h a 3 (— Iuj) 



Inserting + ct 2 — ^)c??7 under the integral sign and making substitutions ai )2 = ^1,2, 
a 3 = rjy, we have A = arj, and the integral in the "radial variable" rj is easily calculated: 



M 3 = 2C F 



9l 



( 47r )d/2 



c<J7 + / dxi dy 5(y + lu'xi — lu) / dr\rfe 



-(XT) 



max(0,w-w') 



where 



uj 



, / /\ w 2 
ur 2 



The integral in y is finite at e — > (and easy to calculate). 

Now we are ready to find the renormalization constant Z + (uj, uj'; //). Re-expressing 



9l 



(47T) d / 

in the matrix element (j6.12j) . we have 

1/2 



M = 7+ Z q Z, 



5(uj - uj') + 2d 



Aire 



x 



6{uj' — uj) uj 



+ 



uj' — UJ UJ' 

6(uj - uj') 
uj — uj' 



+ 5{uj-uj') (6.3i; 

_(l_ log £)^_ w0+o(£) )- 



The quark-field renormalization constants in QCD and HQET th the Feynman gauge are 

Z q = l-C F -±, Z Q = 1 + 2C F -^. (6.32) 



Finally, the renormalization constant "matrix" at one loop is 



Z + (uj, uj'; fi) = 5(uj — uj') + ICp- 



Aire 



x 



6(uj' -uj)uj 6(uj - uj') 



uj — UJ UJ UJ — UJ 

The evolution kernel (J6.9)) has the structure 



+ (-^+ io s^+!|)^-^) 



(6.33) 



T + (uj, uj'; fi) = T(uj, uj'; a s ) + 



UJ 



-T(uj,uj'] a s )log- + 7j(a s ) + l{a s ) 
A 4 



.5(uj-uj'). (6.34) 
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The logarithm log(o;//i) only appears linearly, to all orders of perturbation theory; the 
coefficient of this logarithm is the cusp anomalous dimension [To] : 

T(a s ) = T a s + T ia 2 s + ■■■ T = 4C F (6.35) 

(the two-loop term is also known [15]). The non-logarithmic part is written as jj + 7, 
where the anomalous dimension of the local current j (13. 1|) is 

jj(a s ) = j jo a s + jj!a 2 s H 7 i0 = -3C F (6.36) 

(the two- and three-loop terms are also known |T0J Ej)- It determines the evolution of 
F(/i); the difference 

l(a s ) = 7 a s + 7ia^ H 70 = -2C F (6.37) 

appears in the evolution equation for tp+(uj; [i). The non-5 term is 

r( W y,a.) = r (u,^a. + ... T^M)=( e ^^-, + 6 ^^) ■ (6.38) 

The evolution kernel has dimensionality of 1/energy. 

The evolution equation for the distribution amplitude is 



d log \x 



-r(a s )log- + 7(a s ) 
/i 



(p+(uj; /i) + / r(a;, J\ a s )<p+(uj'; jj)dio' = . (6.39) 



How to solve it? Powers uj n are eigenfunctions of the integral operator in ()6.39j) . by 
dimensionality: 

T(u, J\ a s )J n dJ = f(n, a s )u n , (6.40) 



where 

f (n, a.) = f a s + ■■■ f (n) = 4C F [V>(1 + n) + ^(1 - n) + 2 Ts ] , (6.41) 

from ([6.38]) (here 7^ is the Euler constant). They are not eigenfunctions of the whole 
evolution operator, due to the logarithmic term. We can construct solutions with the 
power depending on //: 



<91og/i \//o7 V/W dlog/i /i 

If the function £ (/i) obeys 



dlog// 

then logcj will cancel in the evolution equation. 
Dividing the definition ([6.42]) by 



<IL r(o.) , (6.42) 



= -2/3(a s ) , /5(a s ) =/5 a s + /3 1 a2 + ... (6.43) 

a log// 
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and integrating, we obtain 



as T(a s ) da s 



i0 2/3(a s ) a s 



2Po 



a s (Y 1 PA 
log — +[ — - — ) {a s - a s0 ) + ■■■ 
a s0 \T fj J 



(6.44) 



where a s o = a s (/io)/(47r). We shall use £ as an independent variable ("time") in the 
evolution equation instead of \i. The function 

f- 



<p(u,0 = - 



MO 



(6.45) 



satisfies the evolution equation 

<9 logy? 



if 



At the leading order, 



lo — + r ( w + £> a ') 

<9£ /i r(a a ) 



(it/ _ _ r(n + g,a 8 ) +7(a g ) 

d£ /i r(a s ) 



(6.46) 



a s = a s0 exp 



and 



_2& 

r ( 

2tt 



log 



// 2vr / 1 



.sO 



2/3 



/i /?o V a 



1 

a s0 



2tt 



— exp — £ 



Vr c 



A) a sO 

g r (n + Q 
r n 



exp —4 



Vr 



1 



Here, from (|6.41|) . 



€ f (n + * = log r(i + n + Qr(i-n) + _ 



r 



T(l - n-£)r(l + n) 



Finally, at the leading order, 



MO 



A \ 2/3n 



r(l -n-£)r(l + n) 

r(i + n + or(i-n) 



exp 



(6.47) 



The distribution amplitude at /^o can be expressed via its moments (|4.20j) : 



+ioo j 

<o;- 1 -">f o) a;-— 



Substituting the solutions ()6.45|) instead of u n gives us the solution of the evolution equa- 
tion (|6.39|) with initial conditions. At the leading order (|6.47p . 



fMi 



exp 



A 



MS 



+oo 



<uj- l - in > { ^ ] u in 



r(l - n - £)r(l + n) dn 
r(i + n + or(i -n) 2^ 



(6.48) 
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Qualitatively, each "bin" in the distribution amplitude at /io near some finite u' pro- 
duces, after evolution to a larger /i, a radiative tail, slowly decreasing (as 1/uj) at u 3> cu', 
due to the evolution kernel T(u,u';a s ) (which behaves as 1/uj at large u, by dimension- 
ality). Therefore, the behaviour of the distribution amplitude at large to is l/tu, up to 
logarithms. The normalization integral 1)4.14)1 of the distribution amplitude logarithmi- 
cally diverges at large ui; its moments (j4.18j) with n > are power-divergent. This means 
that the expression (j4.17|) for the local operators 0+ (n > 0) via + (u) is not valid for the 
renormalized operators. Renormalization of + (u) removes UV divergences in transverse 
momenta; in order to renormalize O ™ , we should also remove longitudinal UV divergences 
(at large to). 

As an illustration let's suppose that at a low fio the distribution amplitude <p+(u)) 
is given by the simple model ()7.24j) . Namely, this initial condition is taken at the scale /zo 
where 



where the quark model is supposed to work, so that the light quark has only momenta of 
order A (no radiative tail). This function is shown by the solid line in Fig. (|T3*j) . where u 
is measured in units of ojq f)7.25j) . Using the leading-order solution ()6.48|) of the evolution 
equation (j6.39J) and supposing A/Aj^g = 1.25, we obtain the distribution amplitude at the 
scale where a s ([i) = 0.5 (dashed line) and 0.3 (dashed-dotted line). We can see that the 
main part of the distribution amplitude (at uj ~ A) becomes lower and the radiative tail 
becomes more prominent when fi grows. 



a s (fj, Q ) = 1 =>- 
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Figure 13: Evolution of _B-meson distribution amplitude <p + (u; n) 
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7 Sum rules 

In order to estimate the quark-antiquark distribution amplitudes ip±(u) of B-meson from 
QCD sum rules, we consider the correlator of the local current j + (}3.2j) (having the quantum 
numbers of the ground-state meson) and the bilocal operator 0±(t) (j4.11j) : 

i<TO ± (t)j+(-x)> = 1± ±^5(x)6(x )fl ± (x°,t). (7.1) 
We are most interested in its Fourier transform 

n ± (*V)= I U ± (x°,t)e^^. (7.2) 

J 27T 

Analytically continuing it from x° > to x° = —ir, we obtain 

n±(r,w) = J p±(e,u)e- £T de = F 2 fi ± (uj)e- AT + U c ± (t,lu), (7.3) 

where p±(e, u) is the spectral density. The correlator contains the contribution of the 
ground-state meson (written explicitly in ()7.3|) ) and of the continuum of excited states. 




Figure 14: Correlator of the local and the bilocal operators 

For sufficiently small r, we can calculate this correlator theoretically. The perturbative 
contribution (Fig. [T4*j) is given in the fixed-point gauge x^A^x) = by the light-quark 
propagator from — x to z: 

H W ( r t ) = - Nc , fl (1) (r, t) = 5 N ° r . (7.4) 

+ v ' ; 27r 2 r(r + lit) 2 ' _ v ' ; 27i 2 r 2 {r + 2it) V ; 

Its Fourier transform is 

n?>(r, u ) = • n ' 1,(T - w) = 4^ e '" T/2 • (7 ' 5) 

Inverting the Laplace transform 1)7.3}) . we find the spectral densities (the integration contour 
should be to the right of the singularity at r = 0): 

p±{e,u>)= / n±(r,o;)e eT — . (7.6) 
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We obtain 

o^sH-D' <£'e.«>=£ (-2) «(*-*)■ <") 

The perturbative spectral density p£ (e, u;) describes the contribution of the on-shell quark- 
antiquark intermediate state with energy e into the correlator (J7.3j) . The on-shell heavy 
antiquark has zero energy; the maximum value of p + of the on-shell light quark with energy 
e is 2e. This explains the ^-functions in ()7.7|) . 

The quark condensate contribution is also important here. It contains the vacuum 
average 

<q(-x)[-x, 0}[0, z]q(z)> . 

This non-collinear quark condensate can be expanded in terms of bilocal condensates ^H] > 
the leading term in this expansion is the bilocal quark condensate 

<q(0)[0,x]q(x)> = <qq> f s (x 2 ) (7.8) 

with x — > x + z. This leading term produces the same contribution into II±: 

fi£V, = -\<qq>fs{{x + zf) . (7.9) 
The expansion of the bilocal quark condensate at small x is 

fs{x 2 ) = l + ^ + ... (7.10) 
lb 

where 

<qG^ v q> = ml<qq> . (7.11) 

In general, it can be written as 

fs(x 2 ) = J fs(v)e vx2 dv, (7.12) 

where fs(u) has the meaning of the virtuality distribution function of vacuum quarks. Its 
moments are expressed via vacuum averages of local operators: 

~ f ~ th ^ 

f s (u)du=l, J f s (u)udu = ^, ... (7.13) 

In terms of this function, the quark-condensate contribution into the Fourier-transformed 
correlators is 

ni 2 V,-) = -^/ s (£) e -^ (7.14) 

In other words, the virtuality distribution function appears directly in the sum rules for 
the 5-meson distribution amplitudes! This is similar to the case of non-diagonal sum rules 
for the pion distribution amplitude [T9*] . 
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The local operator expansion (j7.10j) gives 



fs{") = 8(u)-^S'(v) + :. (7.15) 

Of course, a smooth function can always be expanded in derivatives of (^-function, but such 
an expansion does not tell us much about the shape of this function, unless we sum an 
infinite number of terms. The bilocal quark condensate (J7.8|) at large — x 2 behaves as 

/ s (* 2 )~e~ A ^, (7.16) 

because the Wilson line [0, x] can be considered an HQET heavy-quark propagator, and 
this vacuum average is the correlator of two HQET heavy-light currents, having 5-meson 
as the lowest-energy intermediate state. Of course, the conditions (JZHSJ), (|7.1fij) don't 
determine the shape of the distribution function /s(^) in a unique way. They are satisfied 
by the Bakulev-Mikhailov ansatz 



f s (v) = Ne X p(~-ai?) , (7.17) 

where the parameters a, N are determined by (|7.13j) . 

Now we equate the theoretical result for the correlators to the result obtained from a 
phenomenological model of the spectral densities. They have the contribution ~ 5(e — 
A) of the ground-state meson and that of the continuum of excited states. As usual, 
this contribution is modeled by the perturbative spectral densities ()7.7|) starting from a 
continuum threshold energy e c : 

p±(e, u) = FV±(u;)<$ (e - A) + p£ 1J (e, u)9{e - e c ) . (7.18) 

With this model, the equality of the phenomenological expression for the correlator and 
the theoretical one, 

poo ^oo 

FV±(w)e" Ar + / p { ±\e,uj)e- £T de= / p { ^(e, cu) e - £T de + U ( i\r, u) , 
Jec Jo 

becomes 

FV±Me" Ar = p i l\e,u)e- £T de + U ( i\T,u), (7.19) 
Jo 

where the perturbative spectral density is integrated over the "duality interval" of the 
ground-state meson (from to the continuum threshold e c ). We obtain the sum rules for 
the _B-meson distribution amplitudes 

^ ^ ((- -d T )~ (s) ^ • (7 20) 
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where the functions 

m=0 ' / 

describe the effect of subtracting continuum from the perturbative contribution when the 
spectral density is ~ e n . The perturbative contributions vanish at uj > 2e c , because of the 
properties of the spectral densities (|7.7|) . 

Setting t = in the correlator (J7.1)) . or integrating ()7.2|) in du, we obtain the well-known 
sum rule for F 2 : 

F 2 e~^ = 2^^(s c r) - \<q q >f s {-r 2 ) . (7.22) 
The natural energy scale in this sum rule is 

k= [-— — <qq) ^ 260MeV. (7.23) 



2N, 



c 



With mo/ (4k) = 0.85, one finds the optimal value of the continuum threshold e c /k = 3; 
then a wide plato at l/(kr) G [1.7,2.5] exists, and yields A/k = 1.65. We divide the 
sum rules (|7.20j) by (|7.22j) and take l/{kr) = 2 (in the middle of the plato). The results 
are shown in Fig. where uj is measured in units of k. They are automatically 

normalized (j4.14|) . Of course, the leading-order sum rules cannot tell us at what scale 
these distribution amplitudes are normalized; this scale must be low, \i ~ 1/r. The 
perturbative contributions vanish at w > 2e c . The quark-condensate contribution is the 
same for both ip + (u>) and <p_(a>). Here we used the ansatz ()7.17|) . This contribution gives 
a sharp peak at low u 3 . Details of its shape are unknown, but it cannot be wider because 
of the restriction on the first moment 1)7.13)) . This contribution falls off quickly at larger uj. 
Therefore, the distribution amplitudes at this low \i are only non-zero at u ~ A, in accord 
with the expectations of the quark model. 

If we neglect the quark-condensate contribution (though this is not a good idea) and 
also the effect of continuum subtraction 4 , then the perturbative contributions (j7.5j) suggest 
the following simple model of distribution amplitudes: 

= ^e-^ , <p.(u) = —e-^ . (7.24) 

They are normalized ()4.14j) ; from the first moments (j4.26j) we obtain 

u = h. (7.25) 

These functions (Fig. IT7)) have something in common with the Wandzura-Wilczek ones 
(Fig. EI), but in contrast to them they fall off smoothly at large uj. 

3 For the pion distribution amplitude, the quark-condensate contribution gives enhancements near x — 
and 1. 

4 If the same procedure is applied to the diagonal sum rule for the pion distribution amplitude, it yields 
the asymptotic shape. 
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Figure 15: Distribution amplitude <p+(ou) (solid line), perturbative contribution (dashed 
line), and quark-condensate contribution (dashed-dotted line) 




Figure 16: Distribution amplitude tp-(ou) (solid line), perturbative contribution (dashed 
line), and quark-condensate contribution (dashed-dotted line) 



32 




Figure 17: Model distribution amplitudes: <£+{oj) (solid line) and <f-(oo) (dashed line) 



Radiative corrections to the perturbative spectral density and the dimension-3 quark- 
condensate contribution for </? + were calculated in [Hj. With these corrections, one can 
check that the correlator (|7.2|) satisfies the evolution equation (|6.8|) . The resulting sum 
rules should be used together with the sum rules for F 2 with radiative corrections 
The perturbative spectral density is 



Nc 
8tt 2 



x > 1 : 1 + Cj 



'47T 



+ 2Li 2 



-2 log 2 - - 4(log(x - 1) + l) log - 
— log 2 (x — 1) 



1 — X 



x < 



7 



- (2x + 3)log(x - 1) + 2xloga; + — it 2 + 7 



(7.26) 



x<l: 2C F ^ 
4tt 



2(log(l -x) +x) log - + 21og 2 (l 



x) 



+ (2x — 1) log(l — x) — x 



where x = 2e/u. Now the spectral density does not vanish at e < uj/2, it is just suppressed 
by a s /(47r). Therefore, the perturbative contribution to the sum rule (I7.2()j) for (p+(u>) does 
not vanish at u > 2e c . It produces a radiative tail ~ 1/u with the magnitude of order 
a s /(47r). 

The radiative correction to the dimension-3 quark-condensate contribution cannot be 
used together with a model of bilocal condensate, because radiative corrections to higher- 
dimensional contributions are not known. It is vital to use some model of the bilocal quark 
condensate in the sum rules for the distribution amplitude, because the local operator 
expansion produces contributions 8(u), S'(u)), which don't tell us much about the 
shape of the distribution amplitude. Therefore, the full result for this correction cannot 
be used in the sum rule for ip + {yj). Fortunately, the authors of jH] demonstrated that a 
part of this correction is universal, and exponentialized into the Sudakov factor. We may 
multiply the bilocal quark-condensate contribution by this Sudakov factor. 
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